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Abstract 

We calculate the damping 7 g of collective density oscillations (zero sound) in a one-dimensional 
Fermi gas with dimensionless forward scattering interaction F and quadratic energy dispersion 
k 2 /2m at zero temperature. Using standard many-body perturbation theory, we obtain 7 g from 
the expansion of the inverse irreducible polarization to first order in the effective screened (RPA) 
interaction. For wave-vectors \q\/kF <C F (where kp = mvF is the Fermi wave-vector) we find to 
leading order j q oc \q\ 3 /(vjrm 2 ). On the other hand, for F <C \q\/kp most of the spectral weight 
is carried by the particle-hole continuum, which is distributed over a frequency interval of the 
order of q 2 /m. We also show that zero sound damping leads to a finite maximum proportional to 
| A; — kp\~ 2+2ri of the charge peak in the single-particle spectral function, where n is the anomalous 
dimension. Our prediction agrees with photoemission data for the blue bronze K0.3M0O3. We 
comment on other recent calculations of j q . 

PACS numbers: 71.10.Pm, 71.10.-w 
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I. INTRODUCTION 



The normal metallic state of interacting electrons in one spatial dimension has rather 
exotic properties, which are summarized under the name Luttinger liquid behavior: the 
absence of a discontinuity in the momentum distribution function at the Fermi surface, a 
vanishing density of states at the Fermi energy, and an unusual line-shape of the single- 
particle spectral function, characterized by algebraic singularities and separate peaks for 
spin- and charge excitations^. To discover these features und study them quantitatively, 
it has been extremely useful to have an exactly solvable effective low energy model for in- 
teracting one- dimensional clean metals, the so-called Tomonaga-Luttinger model (TLM)^. 
The exact solubility of the TLM relies on two crucial assumptions: the linearization of the 
energy dispersion for momenta close to the two Fermi momenta ±kp, and the restriction to 
two-body scattering processes involving only momentum transfers small compared with kp 
(forward scattering). The single-particle Green function is then most conveniently obtained 
via bosonization 1 ' 3 . Alternatively, the Ward identity associated with the separate number 
conservation at each Fermi point can be used to derive a closed equation for the single- 
particle Green function^. Within the Ward-identity approach, it is also straightforward to 
show that collective density oscillations can propagate without damping in the TLM, so that 
the random-phase approximation (RPA) for the density- density correlation function H(q, to) 
is exact. As a consequence the dynamic structure factor S(q, u>) = 7r _1 Imn(g, u + iO) of the 
TLM exhibits a sharp ^-function peak. Focusing for simplicity on the spinless TLM with 
forward scattering interactions g 2 = = fo in "g-ology" -notation^, the dynamic structure 
factor is 

StLm(<7,^) = Z q 5(uJ - L0 q ) , (1) 

where Z q = |g|/(27r\/I + F) andw^ = v c \q\. Here F = f /(nvp) is the relevant dimensionless 
interaction, v c = vpy/l + F is the velocity of the collective charge excitations (zero sound, 
abbreviated by ZS from now on), and Vp is the Fermi velocity. 

The simple result (CQ) is a consequence of the approximations inherent in the definition of 
the TLM: the linearization of the energy dispersion and the restriction to forward scattering 
interactions. In more realistic models, we expect that the ZS mode acquires a finite width. 
How does the line-shape of S(q,u) change if we do not linearize the energy dispersion? 
Since the RPA is exact for linear energy dispersion, it is reasonable to use the RPA as a 
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starting point for quadratic dispersion and to try to calculate the corrections to the RPA 
perturbatively. At the first sight it seems that this problem can be solved within the usual 
bosonization approach by treating the effective boson-interactions due to the band curvature 
within conventional many-body perturbation theory for the bosonized problemi^. How- 
ever, for frequencies close to u q this seems not to be possible, and infinite re-summations 
are necessary^. In fact, there are conflicting results for the g-dependence of the damping 7 g 
of the ZS mode of one-dimensional fermions in the literature: while Capurro et alM found 
7 g oc q 3 , Samokhinii, Pustilnik and co-authors^i 3 - 1 ^ and Pereira et at— obtained 7 9 oc q 2 . 
On the other hand, Teber— recently showed that the attenuation rate of an acoustic mode 
embedded in the two-pair continuum of the dynamic structure factor scales as g 3 , which is 
consistent with 7 g oc q 3 . 

Let us briefly discuss the recent works by Pustilnik et al.—^ and by Pereira et al.—, both 
of which found that the damping scales as q . In M . Qit ha, been shown that the width 
7g of the peak of the longitudinal structure factor of the XXZ spin-chain is proportional to 
q 2 . This is not necessarily in conflict with 7 9 oc q 3 for the TLM with band curvature, because 
the XXZ-chain is equivalent to a system of spinless fermions on a lattice. In contrast to the 
forward scattering processes of the TLM, the interaction in this model has also scattering 
processes involving momentum transfers of the order of hp. Although in the Luttinger liquid 
regime of the XXZ-chain these processes are irrelevant in the renormalization group sense, 
their effect on non-universal quantities like 7 g might be essential. A similar argument applies 



also to Ref. 



14j , where the dynamic structure factor of the Calogero- Sutherland model is 



shown to differ from zero only in a finite interval of frequencies of the width proportional 
to q 2 /m. The Fourier transform f q of the interaction in the Calogero-Sutherland model is 
proportional to \q\/m for all momentum transfers q, so that scattering is suppressed for small 
q. On the other hand, the RPA and our strategy of calculating perturbative corrections to 
the RPA can only be expected to be accurate if the interaction f q is strongest for small q. 
We therefore believe that the dynamic structure factor of the Calogero-Sutherland modet^ 
does not represent the generic behavior of the dynamic structure factor of one-dimensional 
Fermi systems with dominant forward scattering, as given by the TLM with curvature. 

As far as the calculation in Ref. [13j is concerned, it is based on an infinite re-summation 
of the apparent leading singularities in the weak coupling expansion, using an analogy with 
the X-ray problem. We believe that this procedure does not properly take into account the 
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FIG. 1: Leading interaction corrections to the irreducible polarization in an expansion in powers 
of the RPA interaction. Solid arrows are non-interacting Green functions and wavy lines denote 
the RPA interaction. 

asymptotic Ward-identity which guarantees the cancellation of all singularities in the limit 
1/m — > 0. Moreover, the renormalization of the real part of the energy of the ZS mode, 
which might be essential to obtain the correct damping^^, has not been properly taken 



into account in Ref. 



13|- 



Within the framework of diagrammatic many-body perturbation theory, the standard 
approach^ to calculate corrections to the RPA is based on the evaluation of the three 
Feynman diagrams shown in Fig. [I], which represent the leading corrections to the irreducible 
polarization in an expansion in powers of the RPA interaction. Surprisingly, these diagrams 
have never been evaluated for the TLM with band curvature, which we shall do in this work. 
It seems that, at least for small interactions, this approach should be sufficient to estimate 
the damping of the ZS mode. We shall further comment on the accuracy of this approach 
in Sec. ED 



II. RPA IN ONE DIMENSION WITH QUADRATIC DISPERSION 

It is instructive to consider first the density- density correlation function for fermions with 
quadratic energy dispersion = k 2 /2m within the RPA, where 

U^ A (Q) = f + Il l (Q). (2) 

For convenience we use the Matsubara formalism and collective labels Q = (iu, q) for wave- 
vector q and bosonic Matsubara frequency iu. To make contact with the usual distinction 
between left-moving and right-moving fermions in the TLM, we write the non-interacting 
polarization as Hq{Q) — Yl a =±^-o(Q)^ where a = + refers to right-moving fermions (with 
velocity Vk = k/m > 0) while a = — denotes left-moving fermions (vk < 0). At finite 
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FIG. 2: RPA dynamic structure factor for quadratic energy dispersion, \q\/kp = 0.1, and different 
values of the dimensionless interaction F. The arrows denote the location of the <5-peak associated 
with the ZS mode, the length of the arrows being proportional to the relative weight W q of the ZS 
peak in the /-sum rule, see Eq. ([5]). 



n o (Q) = -j E Qa ( k ) G o (K)Go(K + Q) , (3) 



temperature T for a system of length L, 

T 
L 

K 

where Gq(K) = [iu — is the free Matsubara Green function and = e a k F +k — £k F = 
avpk + k 2 /2m is the free excitation energy. The label K = (iuj, k) consists of fermionic 
Matsubara frequency iu and momentum k, which is measured relative to akp. The cutoff 
function Q + (k) restricts the range of the fc-integration to the interval [— kp,oo) associated 
with right-moving fermions, while Q~(k) selects the interval (— oo, hp] corresponding to left- 
moving fermions. All degrees of freedom are taken into account in this way. At T = the 
integrations in Eq. (j3J) are easily carried out, 



Tfl 

n° (Q)=a— h, 



avpq — iuj + q 2 /2m 



(4) 



_avpq — iuj — q 2 /2m 

so that dynamic structure factor Srpa^, oj) = 7r _1 Imil RP A(g, to + iQ) can be calculated 
analytically. For any finite value of the interaction Srpa^, consists of two contributions: a 
continuum part due to partice-hole excitations, and a sharp £-peak S^ PA (q, uj) = Z q 5(u — uj q ) 
corresponding to the undamped collective ZS mode, see Fig. El The weight Z q and the 
dispersion u q of the ZS mode are within RPA 

Z q = ^W q , W q = (q/F) 2 /smh 2 (q/F) , (5) 

u q = v F \q\[l + qcoth(q/F) + (q/2F) 2 ] 1 / 2 , (6) 
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FIG. 3: Weight W q of the ZS mode in the /-sum rule in RPA for q/kF = 0.1 (dashed line) as a 
function of F, see Eq. ([5]). 

where q = q/kp. Obviously, the limits of vanishing band curvature (|g| = \q\/mvp —>■ 0) and 
vanishing interaction (F — > 0) do not commute: for \q\ ^> F the weight Z q of the ZS mode 
is exponentially small, so that the particle-hole part of the dynamic structure is dominant. 
For F —>■ the latter approaches a box-function centered at uj = Vp\q\ of width q 2 /m and 
height m/(2ir\q\), see Fig. El In the opposite limit \q\ -C F the ZS peak is dominant. To 
quantify this, note that the dimensionless factor W q defined in Eq. ([5]) can be identified 
with the relative weight of the ZS peak in the /-sum rule J °° duo uj S(q, uj) = Vpq 2 /2n. From 
Fig. [3] it is clear that for \q\ <C F the contribution of the particle-hole continuum to the 
/-sum rule is indeed negligible (actually, it is of order (q/F) 2 and that the crossover 

between the particle-hole regime \q\ ^ F and the ZS regime \q\ F occurs at \q\ ~ F . 

The line-shape of the particle-hole continuum and its position relative to the ZS mode 
are probably not correct within the RPA. In fact, recently Schonhammer— has shown that if 
one uses in the RPA bubbles Hartree-Fock propagators instead of bare ones, then the energy 
of the ZS mode is smaller than the energy of the particle-hole continuum. Moreover, multi- 
pair particle-hole excitations neglected within the RPA will wash out the sharp thresholds 
predicted by the RPA and generate some small spectral weight for all frequencies^. However, 
for sufficiently small q almost the entire spectral weight is carried by the ZS mode, so that 
in this limit we may neglect the particle-hole continuum. 

III. LEADING CORRECTION TO THE RPA 

Because for \q\/k F «C F the spectrum of the density fluctuations is dominated by the 
collective ZS mode, we expect that the exact structure factor for frequencies u close to the 
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exact ZS frequency u q can be approximated by a Lorentzian, 

S{q, u) « ir- l Z qlq [{u - u q f + tJ]" 1 . (7) 

In terms of the irreducible polarization H*{Q) defined via II -1 (Q) = / + n~ x (Q) the inverse 
weight of the ZS peak is given by 

Z~ l = f Q Re%{q,u q ), (8) 

and the damping is 

7q = Imll^g, u q + zO)/RelX(g, w q ), (9) 

where n^(g, u;) = dU*(q, uj)/duj. 

Since 7 g = within RPA, we need to go beyond RPA to estimate 7 ? . This is usually done^ 
by expanding the H*(Q) in powers of the RPA interaction /rpa(Q) = /o[l + fo^oiQ)] 1 - 
Before presenting an explicit calculation, let us give a simple argument for the expected 
g-dependence of j q . Within the functional bosonization approach 7 , U a (Q) can be written as 
a single-particle Green function of a real bosonic quantum field pg representing the density 
fluctuations, 

IT(Q) ex J V[p a }e- s ^p a _ Q p a Q . (10) 

For linear dispersion the effective action S^p"] is quadratic, so that Eq. fflUj) yields the RPA 
result, which is exact for the TLM. Non-linear terms in the energy dispersion renormalize 
the quadratic part of S^p ] and give rise to cubic, quartic, and higher order retarded 
interaction vertices 7 , which all generate corrections to the RPA. We parameterize these 
corrections in terms of an irreducible self-energy E Q (Q). Assuming that Y, a (Q) is analytic 
for small Q (this assumption relies on the cancellation of all singularities in the symmetrized 
closed fermion loops 18 ), and taking into account that Y> a (Q) = Y> a (—Q) (which follows 
from the invariance of the right-hand side of Eq. fflQj) under Q — > —Q), we conclude that 
E a (g, uj q + iO) = c + C2(q/kp) 2 + 0(g 4 ), with real c and complex c 2 . But the relation 
between the density fields pq and the Tomonaga-Luttinger bosons whose damping 7 q we are 
seeking involves an extra factor of \q\ 1 ^ 2 (see for example p. 57 of Ref. JtJ), so that we expect 
7q oc |g|ImS a (g, uj q + iO) oc \q\ 3 . 

We now confirm this result by explicitly calculating T, a (Q) to first order in the RPA 
interaction. The diagrams contributing to Tl*(Q) to this order are shown in Fig. [TJ Writing 
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again n„(Q) = £ a n?(Q) we have n?(Q) « n?(Q)+n? (Q), with n?(Q) = n« (Q)+n« (Q). 
The sum of the self-energy corrections shown in Fig. [1] (a,b) is 



L 2 

JsT.Q' 

xG«(^ + Q') [Gfttf + Q) + - Q)\ , 
and the vertex correction in Fig. [1] (c) is 



(11) 



nf B (g) 



L 2 



K,Q' 



xG«(K + Q')G«(K + Q + Q') 



(12) 



To regularize some of the q'- integrations, we introduce a momentum transfer cutoff q c ^kp 
which restricts the integration range to < q c . The imaginary part of these expressions is 
not sensitive to q c . To evaluate Eqs. ( TlTlfT2l) in the limit T — > it is convenient to write^ 

f 2u'S RPA (q,u') 



poo 

/bpa(9, «*>) = /o - /o / da;' 

Jo 



(13) 



a/ 2 + cu 2 ' 

The frequency integrations in Eqs. fllip and 0121) can then be carried out exactly. The result 



is 



c-2 
_0_ 

L 2 



L 2 

k,q' 

OO f 



+ 



sgn(#)0(-ffff +q ) 
[4° - £ + , + M 2 



^0 a (A;) y dcj' 

k,q' o 

Srpa(9W) 



8 - H% +q + 



a ■> 



[4-4% + M 2 



141 + 14V + ~ + iffi^j + W + sgn(4 a +g , 

-> (-<?, -iu;)] , 



17 V SfcSfc+«j' 



-£j^)sgn(#)0(-tf£j^) 



Srpa(<?W) 



X 



(14) 



(15) 
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For linearized dispersion, av F k, we have verified that II" s (g, zo;) + IT^g, ia;) = 0, in 

agreement with the closed loop theorem^ 7 . With quadratic dispersion = ctt^fe + k 2 /2m 
this cancellation is not perfect. Assuming that the band curvature 1/m is small, we may 
approximate Srpa(9S ~ Srpa^, <^)|i/m=o = <Stlm(<7) in Eqs. f[T4"l) and (fT5|) . because 
the m-dependence of 5 , rpa(<?, u;) is irrelevant for the cancellation between self-energy and 
vertex corrections in the limit 1/m — > 0. One should keep in mind, however, that this is 
only a good approximation as long as the dynamic structure factor is dominated by the ZS 
mode, which is the case for \q\/kp <C F. We shall from now on focus on this regime. With 
<Stlm(<Z, oj) given in Eq. (CQ), the frequency integration is trivial. However, the resulting 
expression is ill-defined, because IrnII"(g, u + iO) contains singular terms proportional to 
5(uj — The fact that the direct expansion of II* (g, oj) in powers of the interaction 

generates unphysical singularities has been noticed long time ago^. These singularities can 
be avoided by expanding the inverse polarization [II" (Q)] -1 in powers of the interaction, 
[II^(Q)]- 1 = [U^Q)}- 1 - where to first order in the RPA interaction T%(Q) = 

[n^(Q)]" 2 n°(Q), with n°(Q) = nf s (g) + n^(Q) given in Eqs. PJ and The point is 
that according to Eq. ( |T0l) U a (Q) can be viewed as a bosonic single-particle Green function, 
which should never be calculated directly, because its expansion in powers of the interaction 
usually contains unphysical singularities associated with the free Green function. On the 
other hand, the expansion of the inverse Green function (i.e., the self-energy) is expected to 
be regular. Ignoring the renormalization of the real part of the ZS dispersion, we approximate 

7 g ^ImS^(g,^ + z0)[n^(g,^)] 2 /n[ ) (g,^) , (16) 

where a q = sgn(g). For simplicity we consider only the leading behavior of 7 g in an expansion 
in powers of q and 1/m. The integrations over the momenta k and q' in Eqs. (1141) and ( |T5l) 
can then be carried out exactly. For \q\/kp min{F, 1} we obtain the expected result 

_ 7T \(f_ 

lq ~ 8 v^TF[l + VTTF} A v F m 2 ' { } 

An intensity plot of the corresponding structure factor (J7J) in Lorentzian approximation is 
shown in Fig. HI 

Our result agrees qualitatively with Ref. [h]], which found 7 9 oc F\q\ 3 /vprn 2 , but we 



disagree with Refs. 1 ill 121 ]. who obtained 7 g oc q 2 /m. Possibly the discrepancy with Eq. ([T 



is related to the non-commutativity of the limits q — > and F — > 0, which is obvious from 




FIG. 4: Dynamic structure factor in Lorentzian approximation for F = 1, see Eqs. (I5|6|7p and 

(HZD- 




-1.75-1.5-1.25 -1 -0.75-0.5 

H\q\/k F ) 

FIG. 5: Maximum intensity I max (q) of the charge peak in the photoemission spectrum of the quasi 
one-dimensional metal K0.3M0O3. The points are data from Fig. 8 of Ref.— , assuming that q = 
corresponds to an emission angle 9 = 7°. The lines are fits to I m ax(<z) oc |g| _2+2r '. Our theory is 
valid only for \q\ <C kp, which might explain the discrepancy between positive and negative q for 
larger \q\. 

Eqs. (jSJ) and (jSJ). Our result 7 g oc q 3 is valid for \q\/k F <C F, which for fixed \q\ requires a 
minimum strength F of the interaction, but for fixed F > is always satisfied for sufficiently 
small q. In the opposite limit F \q\/kp, where most of the spectral weight is carried by 
single particle-hole excitations, our approach also yields 7 g oc q 2 /m. The g 3 -scaling of 7 g 
is consistent with a recent result by Teber 16 , who showed that the attenuation rate of a 
coherent acoustic mode embedded in the two-pair excitation continuum scales as q 3 . 

S(q, uj) can be measured using X-ray scattering. Unfortunately, for small q the intensity 
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is very low, so that an experimental verification of Fig. H] seems to be difficult. However, 
the ZS damping has a dramatic effect on the shape of the single-particle spectral function 
A(kp + q, u) for u pa u q , which can be measured via photoemission. For the TLM (including 
now the spin degree of freedom) one finds for uj — > u q = v c \q\ an algebraic singularity^, 
A{kp + q,ui) oc |g(w-w g )| (lJ-1)/2 , which is a consequence of the undamped ZS propagation. 
Here r\ is the anomalous dimension. The damping of the ZS mode washes out this singularity. 
The sensitivity of the spectral line-shape for u ~ uj q to non-universal perturbations neglected 
in the TLM has previously been noticed by Meden 2 ^. In fact, from the expression for 
the single-particle Green function derived via functional bosonization^ we estimate that 
the photoemission intensity exhibits for \uo — uo q \ <C 7 9 a finite 2 ^ maximum, I max (q) oc 

In Fig. we compare this prediction with photoemission data^ 4 - for 
the "blue bronze" K0.3M0O3. The resulting value of r] is consistent with previous estimates 24 
between 0.7 and 0.9 for this material. However, the dataM were taken at T = 300A", so that 
deviations from our T = theory are expected 2 ^. Refined photoemission data probing the 
T = regime of a Luttinger liquid would be useful. 



IV. SUMMARY AND CONCLUSIONS 



In summary, using a standard diagrammatic approach we have calculated the damping 
7 9 of the collective charge mode (zero sound) in a Luttinger liquid due to the non-linearity 
in the energy dispersion. Our result j q oc \q\ 3 /v pm 2 suggests that ZS is indeed a well defined 
elementary excitation in a Luttinger liquid. We have pointed out that the photoemission 
line-shape close to the charge peak is very sensitive to the ZS damping and have made a 
prediction for the height of the charge peak which agrees with experiment^ 4 -. ZS damping is 
also crucial to understand Coulomb drag experiments in quantum wires 12 . 

We emphasize that our result is based on the expansion of the inverse irreducible po- 
larization to first order in the RPA interaction, using standard many-body perturbation 
theory. For Fermi systems in three dimensions this approach has been quite successful 1 ^, 
and we have shown here that also in one dimension the perturbative calculation of the ZS 
damping does not give rise to any singularities, at least to first order in an expansion in 



powers of the RPA interaction. However, according to Ref. 



13| higher order corrections in 



the bare interaction (not contained in our expansion to first order in the RPA interaction) 
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generate new singularities which, when properly re-summed, lead to a much larger damping, 
7g oc q 2 /m. We suspect that the re-summation procedure proposed in Ref. [3] is unreliable, 
because it does not properly take into account the closed loop theorem^. We are currently 
investigating this problem using a functional renormalization group approach^ 1 ^ with mo- 
mentum transfer cutoff, where corrections to the RPA can be calculated systematically even 
in the presence of infrared singularities. 

Finally, let us emphasize that, in contrast to expansions using conventional 
bosonization^&ii, our approach is not based on a direct expansion of the dynamic structure 
factor in powers of 1/m. This is obvious from the fact that in the non-interacting limit we 
recover the exact dynamic structure factor of the free Fermi gas with quadratic energy dis- 
persion. Our approach can be formally justified with the help of functional bosonization^^, 
where the vertices of the effective bosonized theory are given by symmetrized closed fermion 
loops and the effective small parameter which controls a perturbative expansion is propor- 
tional to the combination Fq c /kp, where the range q c of the interaction in momentum space 
must be small compared with kp — mvp. 
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